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Abstract 

We study the nonlinear cr-model in an external magnetic field applied on curved surfaces with 
rotational symmetry. The Euler-Lagrange equations derived from the Hamiltonian yield the dou- 
ble sine-Gordon equation (DSG) provided the magnetic field is tuned with the curvature of the 
surface. A 2tt skyrmion appears like a solution for this model and surface deformations are pre- 
dicted at the sector where the spins point in the opposite direction to the magnetic field. We also 
study some specific examples by applying the model on three rotationally symmetric surfaces: 
the cylinder, the catenoid and the hyperboloid. 
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The interplay between geometry and physical properties of condensed matter systems has 
demanded large attention in last decades. In part, this is associated to the growing capacity 
for fabricating and manipulating nanoscale devices with different geometries. For instance, re- 
cent developments in nanotechnology made it possible to produce quasi two-dimensional exotic 
shapes like the Mobius stripe JH, torus J2l and asymetric nanorings 10. On the other hand, there 
is a growing theoretical interest in this subject due the fact that particle-like excitations, like vor- 
tices and skyrmions, for example, interact not only with each other, but also with the curvature 
of the substract 10]. Vortices and skyrmions appear in several condensed matter systems, e.g. 
superconductors, BEC, nematic liquid crystals, ferro and antiferromagnets 05l|6|]. In the case of 
ferromagnetic materials, cylindrical nanomagnets with a vortex as the magnetization groundstate 
have been considered as candidates to be used in logic memory, data storage, highly sensitive 
sensors [7] and cancer therapy [8]. The energy of these topological excitations depends on the 



curvature of the ferromagnetic nanoparticle Q9l |10| and their dynamical properties are affected 
by the interaction with curve defects appearing during the fabrication of the nanomagnets ll ill . 

In this context, an important issue is the control and manipulation of the shape and physical 
properties of nanoparticles. Indeed, several works have already devoted attention on this issue, 
e.g., theoretical works have shown that the application of a constant external magnetic field in 
a circular elastic cylinder surface coated by magnetic material must cause a deformation at the 
region where the spins are pointing in the opposite direction of the magnetic field 11211 . At the 
same time, the concept of magnetoelastic metamaterials, which respond by being compressed by 
electromagnetic forces in the presence of an external magnetic or electric field, leads to a new 
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generation of metamaterials that can be useful for the design of artificial media [13]. In addition, 
it was shown that the curvature of graphene bubbles can be controlled by applying an electric field 



1 14j] and still, by combining the curvature effects with magnetic fields, the molecular alignment 
of flux lines of the nematic director of nematic liquid crystals can be reoriented or switched 
between two stable configurations [15]. 

In this letter, we show that the homogeneous DSG can be obtained for an arbitrary mag- 
netically coated surface with rotational symmetry in the presence of an external magnetic field 
provided the field is tuned with the curvature of the surface. The obtained solution consists 
in a 2n skyrmion, which induces a deformation in the rotationally symmetric surfaces due the 
presence of the magnetic field. Besides, we apply the model to three specific cases: cylinder, 
catenoid and hyperboloid, in order to obtain the characteristic lengths of the skyrmion for these 
particular surfaces. To our knowledge, this is the first time that the tuning between an exter- 
nal magnetic/electric field and the geometry of the substract yields a controllable system with 
well defined theoretical solutions. The coupling between the magnetic field and the geometry 
of condensed matter systems is a surprising result and we believe that this issue could guide 
future studies in the control and manipulation of the shape and physical properties of magnetoe- 
lastic coated surfaces, superfluid helium, nematic liquid crystals, graphene bubbles, topological 
insulators and so forth. 

The total energy for a deformable, magnetoelastically coupled manifold is given by E — 
E mag +E e i+E m . e \, where E, mg , E e \ and E m . e \ are the magnetic, elastic and magnetoelastic contribu- 
tions to the total energy, respectively. The magnetic contribution is composed by the exchange, 
magnetostatic and Zeeman terms. In this work, we will focus our attention mainly on the ex- 
change and Zeeman contributions to describe the magnetic properties of curved surfaces with 
rotational symmetry in the presence of an external magnetic field. In this case, the energy is well 
represented by the nonlinear cr-model on a surface in an external axial magnetic field: 



// = jj (VmfdS-gfi jj m -BdS, 



(1) 



where m is the magnetization unit vector, dS is the surface element, B is the applied magnetic 
field, jj is the magnetic moment, and g is the g factor of the electrons in the magnetic material. 

This model was previously considered for circular cylinder surface and the well known DSG 
lfl6ll was derived from the Euler-Lagrange equations lll2"l . In that case, the authors obtained a 2n 
skyrmion like solution and showed that a geometrical frustration appears due to a second length 
scale introduced by the magnetic field. On the other hand, unlike the cylinder case, when the 
above model was considered on an arbitrary curved surface, it yielded the inhomogeneous DSG, 
whose solution could be gotten only numerically |l7ll. Here, in order to obtain the homogeneous 
DSG system, we will consider that the magnetic field is in the z-axis direction and its strength is 
a function of p, that is, B = B(p)z, where p = p(z) is the radial distance from one point on the 
surface to the z axis. 

To proceed with our analysis, it will be useful to rewrite Eq. (Q]) on a general geometry with 
metric tensor gjj. In this case, the Hamiltonian is given by: 



H= ff ^^k^ 5 p ) 2+ ^ sin20 M 2+ ^x i - cose) ) 



dpdcp, (2) 



where g pp and g^ are the contravariant metric elements. The magnetization unit vector is 
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parametrized by m = (sin cos <t>, sin © sin <1>, cos 0) and the rotationally symmetric curved sur- 
faces are parametrized in cylindrical coordinates system (p, <p, z). 

The Euler-Lagrange equations for the Hamiltonian (|2} give us the following equations: 




and 




sin20+— B'(p)sin0 



= 0, 



(3) 



(4) 



where B'(p) = gfiB(p). 

Once we are considering surfaces with rotational symmetry, the parametric equations associ- 
ated to these can be given by r = [p cos (p, p sin <f>, zip)), where p is the radius of the surface at 
height z, and <p accounts for the azimuthal angle. In this case, we have that the covariant metric 
elements are given by: 



= J_ = 2 



and 



8pp = 



\dpj 



what implies that d<p ^ -y^l - 0. Thus, the Eq.© leads to: 

d 2 ,® = -» O = 4> + 0o, 



(5) 



(6) 



where 0o is a constant of integration that does not influences the energy calculations. Now, we 
can rewrite Eq.© as: 



d)® = 



sin 20 



+ g<p<pB'(p) sin 0, 



(7) 



where di; = y/gW/gPPdp. 

From the Eq.©, there are two situations to be considered: firstly, when we consider B(p) = 0, 
it yields the sine-Gordon equation and a skyrmion, which satisfies the self-dual equations, is 
obtained II 1 811 - T his case was previously studied for many different curved geometries, e.g., 
cylinder B18lll9ll . sphere |120|1, pseudosphere [21], torus 112211 . cone 12311 . catenoid, hyperboloid 
iflvL fToh . It is important to note that, usually, the skyrmion presents a characteristic length scale 
(CLS) that depends on the length of the underlying manifold, which, in general, appears in front 
of the sin(20) term. However, in our model, this characteristic length is given by ( sfg^ /g pp ), 
which is embedded in the £ parameter. Thus, the skyrmion would have its CLS rescheduled to 
one, that is, if we take the cylinder case, the CLS of the skyrmion is equal to the radius p B of 
the cylinder. However, when working on an infinite cylinder, the change of variable z — > z/p 
eliminates any dependence on p lfl8ll . 

The second case to be analysed is given when we consider an arbitrary magnetic field B(p) + 
0. In this case, Eq.((7]) gives the inhomogeneous DSG and numerical solutions are demanded 
lfl7ll . However, note that if the external field is tuned with the surface curvature in the form 
B'(p) = g^BL where B' = gpB(p ) and p is the surface radius at z = plane, Eq.(Q yields the 
homogeneous DSG: 

d]&= ^^+B o sln0, (8) 



whose solution can be given by: 



©(£) = 2 tan" 



£ sinh | 1 



(9) 



where p 2 B = l/B' Q and £ = p a /(l + p 2 ) 1 ^ 2 - Then, one can see that the homogeneous DSG must 
be obtained provided that the magnetic field is tuned to the surface geometry, depending on the 
radius p. 

The above equation represents a 2n skyrmion, which is a topological excitation belonging 
to the second class of the second homotopy group, whose CLS is given by Note that the 
increasing of B' Q leads to the decreasing of the CLS of the skyrmion causing it to remain confined 
in smaller regions of the surface. It is easy to note that £ — > when B' — > oo. See also that if 



we do B' 



in Eq.©, we do not get the n skyrmion, as it happens when we consider Eq.(|8} 
with B' = 0. This is explained because these two solutions belong to two different homotopy 
classes and cannot be transformed one to another by continuous transformation or by limiting 
process B' Q — > oo. By introducing a CLS in the system, the magnetic field induces a geometrical 
frustration, and the skyrmion, which had its CLS rescheduled to one, must choose a new CLS, 
given by that is smaller than the length p B , introduced by the magnetic field. 

In order to calculate the energy of the 2n skyrmion with the imposed conditions, we can rewrite 
the Hamiltonian (0 as: 



H' = Hi + H 2 = 2n 



I [(<9 f ©) 2 + sin 2 ©] d£ + I — ( 1 - cos 0) 

tj — CO *J —oo P B 



d£, 



(10) 



where Hi is the part of the Hamiltonian that does not have dependence on the magnetic field and 
H<i is the Zeeman term. In this case, the energy of the skyrmion is evaluated to give: 



£ r = 8tt 



1 \ 1/2 1 
1 + — + — sinh 1 p B 

P B ) P B 



(11) 



that is larger than the minimum energy for the homotopy class with winding number Q — 2, 
that is to say, E 2nS = 2>n. Only in the limit p B — > oo, we get E s — > 87r and ( — > 1. Once 
p B = ^1/B'(p ), the energy of the skyrmion is associated with the strength of the magnetic field 
at the z — plane. In order to release the geometrical frustration introduced by the magnetic 
field, an elastic surface will deform to decrease the radius in the region where the skyrmion is 
centered. This will minimize the energy of the two terms that does not have dependence on the 
magnetic field in the Hamiltonian ( TTUb . 

Now, we will focus on the last term in the Hamiltonian H', which involve the external magnetic 
field: 



f 00 1 
H 2 = 2n - 2 {l- 

For the solution ©, Hi will take the form: 



cos ®)d£. 



H 2 = 4n 



r° 1 

\J — oo P B 



1 



1 +(1 -C 2 ) sinh 



2 i 
I 



<J —i 



(12) 



(13) 



This part of the Hamiltonian will decrease in energy if we locally decrease £, but keep it constant 
in Eq.©. Therefore, for a magnetic field tuned with the surface, if we keep the surface cross 
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8tt 



1 + — (l+21n2pj 



(14) 



section at £ + oo, the surface will deform in the region of the skyrmion. In the case of small 
magnetic field, we get: 

1 

and for large magnetic field, we get E s = 8^/p s . 

The highest contribution for the magnetic energy density due the interaction with the external 
magnetic field comes from the sector where the spins point in the opposite direction to the field. 
In this way, the 2n skyrmion would like to collapse and eliminate the region where the spins are 
opposite to the magnetic field. However, to align the spins, the system would fall in a n skyrmion 
sector of the second homotopy group, which has infinite energy due the interaction of the spins 
at +00 with the magnetic field. Thus, there is a hard-core repulsion between the two n skyrmions 
in the system by virtue of the curvature. Due the introduction of a geometric frustration in the 
problem by the magnetic field, the energy in Eq. dTTb is grater than 87r. Then, the surface will 
deform to decrease its energy till it have that one belonging to the second homotopy class 01211 . 

In the absence of an external field (B = 0), the self dual equation for the Eq.© is <9^0 - 
+ sin 0. For the DSG, we can write: 



d f = + sin0 



i + 4sin 2 (0/2) 



p 2 sin 



(15) 



The function in square brackets is greater than that one for all values of 0(£) and the self dual 
equation is not satisfied. Therefore, the surface would deform in order to satisfy the self dual 
equation and to lower the energy. Note that all solutions of the Eq.(fT5l) satisfy the Euler-Lagrange 
equation, but not vice versa. The 2n skyrmion lattice of the double sine-Gordon equation in two 
regimes, p B < 1 and p B > 1, are similar to that found on the cylinder II 1211 . 

Now, we will apply our generic model in some surfaces with rotational symmetry in order to 
find the skyrmion CLS associated to each one. Firstly, we will investigate the above model on 
the cylinder surface, which was previously studied by Saxena et al II 1 211 . In this case, we have 
that zip) = z is a constant. This parametrization leads to g c J. = p 2 = r 2 and g C pp = 1, where 
r is the radius of the cylinder. Thus the magnetic field to be applied on the surface would be 
B'{p) = (l/r 2 )B' Q . Since r is constant, the external magnetic field is also a constant. Finally, we 
have that f . = In r. As expected, the result obtained here would be the same as that given in the 



Eq.(4) of Ref. 1 12] if we had considered z instead p as the variable of the problem. 

To our knowledge, unlike the cylinder, the homogeneous DSG for Heisenberg spins in an 
external field was not found in any other surface and attempts to find solutions to the model 
described by the equation (Q} led to the inhomogeneous DSG, which must be solved numerically 
lfl7ll . Here, we will focus our attention on two surfaces with negative varying Gaussian curvature, 
the catenoid and the hyperboloid. The catenoid is a minimal surface, as its mean curvature is 
zero everywhere. For this surface, we have that p(z) = r cosh(z/r), thus, this surface can be 
parametrized by (pcos 0,psin0, rcosh _1 (p/r)), where r is the radius of the catenoid at the z — 
plane. This parametrization leads to: 

Then, the magnetic field to be applied on the surface, in order to obtain the homogeneous DSG, 
will be B'(p) m = ( l/p 2 )B' . Once p is not a constant for the catenoid, unlike the result obtained 
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for the cylinder, the field must not be constant. Indeed, for p — > oo, the field tends to zero and 
its maximum value (B'(p = r) cal = (l/r 2 )B' Q ) must occur at p — r. Finally, the catenoid in an 
external varying magnetic field admits a 2n skyrmion solution with 



(17) 



The hyperboloid has a shape similar to the catenoid, however, while the catenoid has mean 
curvature null everywhere, the hyperboloid has variable mean and Gaussian curvatures. A hy- 
perboloid can be parametrized by (pcos0,psin</>, (b/r) yjp 2 - r 2 ), where again r is the radius of 
the surface at z — and b is a multiplicative parameter that accounts for the height of the surface. 
This parametrization yields: 

h p 2 (b 2 + r 2 )-r 4 
fi>=t and g %= )—±- . (18) 

From now on, we will use the biharmonic coordinate system (BC), in which b — r, to describe a 
particular kind of hyperboloid ll24ll . so called polar hyperboloid. In this case, the metric element 
is simplified to: 

Spp YP = (19) 

and, to obtain the homogeneous DSG, the magnetic field must be given by B'(p) hyr = (l/p 2 )B[ y 
Thus, as well as the catenoid case, the field must be p dependent, having its maximum value, 
given by B'(p = r) hyp = (l/r 2 )B' , at the plane z — 0. The magnetic field applied on the hyper- 
boloid described by BC varies with p = Vr 2 + z 2 . Then, B'(p) hyp — > when z(p) — > ±°°. Finally, 
the polar hyperboloid also admits the solution given in the equation (0, however, the £ parameter 
for this surface is given by a large and tedious expression in such a way that it will be omitted 
here. Despite the fact that linip^co B(p) = 0, it is important to note that the total magnetic flux 
through these infinite surfaces diverges as In p at infinity for the catenoid and hyperboloid, but 
remains constant and depending on r for the cylinder. 

In conclusion, we have shown that the Euler-Lagrange equations derived from the continuum 
approximation for classical Heisenberg spins on a surface with rotational symmetry in an external 
magnetic field in the z axis direction is the homogeneous DSG provided the field is tuned with 
the curvature of the surface. We have found a single 2n skyrmion-like solution to this model for 
an arbitrary surface and surface deformations at the sector where the spins point in the opposite 
direction of the magnetic field were predicted for an elastic surface. The model was applied 
to three specific cases: the cylinder, the catenoid and the hyperboloid, which have different 
characteristic length scales and curvatures. Each considered surface admits a skyrmion solution, 
but the magnetic field which gives the homogeneous DSG presents different behavior to each 
one, once it is a function of p(z). While the field on the cylinder is constant, that one which 
leads to the DSG on the catenoid and the hyperboloid is inversely proportional to p(z). Since the 
characteristic length of the skyrmion depends on the strength of the field, the magnitude of the 
deformation can be controlled by the magnitude of the applied magnetic field. 

These results may play an important role for experimental studies where the magnetoelastic 
effects predicted here could be observed. Besides, theoretical works including surface tension 
could clarify the deformation dynamics of magnetic coated membranes and/or magnetoelastic 



metamaterials 01311 . Finally, we believe that this issue may be relevant for future studies on 
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manipulation and control of the shape and physical properties of curved graphene sheets and 
topological insulators. 
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